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Stability of the curvature perturbation in dark sectors' mutual interacting models 
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Abstract 

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove 
that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy 
is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the 
curvature perturbations; however if the coupling is proportional to the energy density of dark energy, there is room for the 
stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no 
matter whether the coupling is proportional to the one or the other energy density. 
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We are convinced by the fact that our universe is undergoing an accelerated expansion driven by the so called dark 
energy (DE) . The leading interpretation of such a DE is the cosmological constant with equation of state (EoS) w = —1. 
Although the cosmological constant is consistent with the observational data, it presents apparently unsurmountable 
problems from the theoretical point of view, which not only requires a severe fine tuning of 120 digits to attain the 
actual value of the cosmological constant, but also leads to the coincidence problem, namely why the DE and the 
Dark Matter (DM) are comparable in size exactly today [1]. 

Dark Energy contributes a significant fraction of the content of the universe. It is thus natural to consider its 
interaction with the remaining fields of the Standard Model in the framework of standard field theory. The possibility 
that DE and DM can interact has been widely discussed recently [2]-[27]. It has been shown that certain types of 
coupling between DE and DM can lead to a late time attractor solution for the ratio of DM and DE densities [()] and 
provide a mechanism to alleviate the coincidence problem [2, 4]. It has been argued that an appropriate interaction 
between DE and DM can influence the perturbation dynamics and affect the lowest multipoles of the CMB spectrum 
[9, 11]. Arguments using galaxies structure formation suggested that the strength of the couphng could be as large as 
the QED fine structure constant [9, 12]. More recently, it was shown that such an interaction could be inferred from 
the expansion history of the universe, as manifested in, e.g., the supernova data together with CMB and large-scale 
structure information[l(), 23, 24]. In addition, it was suggested that the dynamical equilibrium of collapsed structures 
can be affected by the coupling between DE and DM [13]. The basic idea is that the virial theorem is distorted by 
the non-conservation of mass caused by the coupling[22]. Thcrmodynamical attempts to understand the interaction 
between DE and DM has also been proposed [18]. 

Recently there has been some concern about the stability of the perturbations under DE and DM interaction [27], 
which could represent a sharp blade in the heart of such interacting models. In the original analysis the authors 
considered that the energy exchange between DE and DM is proportional to the energy densities of DM and total 
dark sectors. For the constant DE EOS w > —1, it was found that the instability arises regardless of how weak the 
coupling is. In this work we are going to reexamine the stability of the curvature perturbation when dark sectors are 
mutually interacting. We will concentrate on the interaction between dark sectors in a linear combination of energy 
densities of DE and DM, which is a more general phenomenological form in describing the interaction [23, 29]. We 
will restrict our investigation to constant EOS including w > —1 and w < —I cases. 

We consider a two-component system with each energy-momentum tensor satisfying 

v.r^';'; = g;^,, (i) 

where Q'rx) denotes the interaction between different components and A denotes either the DM or the DE sector. This 
equation can be projected on the time or on the space direction of the comoving observer. Using the four velocity V^, 
it can be contracted into 

K v^T^';-; = -px-9{p,+px) = KO^A) > (2) 

which is the projection on the time direction of the comoving observer. Above, p\ ~ V^V^ px and 6 ~ V^Vi, is the 
volume expansion rate. In order to get the projection along the space direction, we can use h^ ~ S'^ + V^Vu on (1) 
and take the contraction 

^: v^rf^-; = [px + vx)A-' +^'^ V>A = KQlx) ■ (3) 

where A'^ ~ V^Vp^V^ is the acceleration. For the homogeneous and isotropic universe, it requires '■^■'V'^pa = 0. 
Besides, for the DM particle, its world line is the geodesic, A'^ = 0. Thus the spacial part of Q'^^. vanishes, which 
means that in the background there is no momentum transfer between dark sectors [27]. For the whole system the 
energy momentum conservation still holds, satisfying EaV^T/^"^ = 0, thus requiring Q%^ = —Q%m- 
We choose the perturbed space-time 

ds^ = a^[-(l + 2i;)dT^ + 2diBdTdx' + (1 + 2(t>)5^jdx'dx^ + D.jEdx'dx^], (4) 

where 

A, = (d.d, - i<5.,v2). (5) 



The perturbed energy-momentum tensor reads 



sv.Ti;^^ 
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regardless of the anisotropic stress, while (5Qk •, is a new perturbation variable. Considering that the intrinsic momen- 
tum transfer can produce acoustics in the DM fluid as well as pressure which may resist the attraction of gravity and 
hinder the growth of gravity fluctuations during tightly coupled photon baryon period, in our study we shall neglect 
the intrinsic momentum transfer by setting SQl-^-. = 0. This is a choice of interaction and the results should not 
heavily depend on such assumption. Our aim is to provide examples of both stability and instability in perturbations. 
We construct gauge- invariant quantities by employing Bardeen's potentials, gauge-invariant density contrast and 
velocity 
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Choosing a particular gauge, the Longitudinal gauge, by defining i? = 0, i? = 0, one can find 'I' = -0, <J> = (/) [2n] 
For the interacting model we use the perturbed pressure [27] 



5pd = ClSdPd + [Cl - Cl) 



3H(1 + w)VdPd 2n0^d 
k "^"^T 



and the interaction as a linear combination of the energy densities of dark sectors, 

a'^Q'ii = 3H{Xip„i + \2Pd) 
a^Q" = -37Y(AiPm -I- \2Pd) 



(7) 



(8) 



(9) 



where Ai and A2 are small positive dimensionless constants. The generality in the choice of the couplings relies on the 
generality of the models. In case we had a Lagrangian formulation the coupling should be fixed. Lacking a Lagrangian 
we are free to choose our model. We are going to show that for some choice of couplings we may achieve stability in 
curvature perturbations. Choosing a positive sign for the interaction the direction of the energy transfer goes from 
DE to DM, which is required to alleviate the coincidence problem [15] and avoid some unphysical problems such as 
negative DE density etc[23, 27]. In [27] it was argued that it is more natural to assume that the interaction between 
dark sectors depends on purely local quantities. Considering the symmetries of the Friedmann-Robertson- Walker 
metric, we note that the interaction can vary only in time, rather than from point to point. The only time parameter 
in question is the age. Thus, the factor Ti. appears in our interaction, implying that the interaction depends on the 
cosmic time through the global expansion rate. 

By taking Fourier transformation of eq( 6), we get perturbation equations 



D'^^ = -kU,n + m^{\i + \2/r)~i{\i + X2/r)^' + m\2{Dd- Dra)/r 

U'„, = -nUra + k^-m{\l+X2lr)Um , 



D'^ = -:mcl{Dd-[i{\ir + \2) + i{i + w)]^--m{cl-cl) 



mUd 



Ud 



k ''^''{l+w)pdk 



U'd 



-iT-Lw [3(Air + A2) + 3(1 + w)] $ + -iHwDd + 3w'$ + 3(Air + As)*' - kUd - 6*H(Air + A2) 

+mXir{Dd - D„,) 

-nil - 3w)Ud + kCl {Dd - 3[(Air + A2) + (1 + u;)]$} 



<Cl~Cl)a^Q'd 



Ud 



(1 + w)pd 



i{Cl ~ CDnUd + (1 + w)fc* + 3H(Air + \2)Ud. 



(10) 
(11) 



(12) 
(13) 



where r ~ Pm/pd, U = {1 + w)V. We have taken SH = by assuming the expansion rate in the interaction to be the 
global expansion rate. This is a matter of choice. 

In the above, C^ = w < 0. However, it is not clear what expression should we have for Cg. In [27] it has been 
argued in favor of Cg = 1. This is correct for the scalar field, but it is not obvious for other cases, especially for 
a fluid with a constant equation of state. The most dangerous possibility, as far as instabilities are concerned, is 
C^ ~ 1 ^ C^ ~ w < since the first term in the second line of eq( 13) can lead to blow up when w close to —1. In 
spite of such a danger we arc considering such a case here. Assuming Cg = 1, C^ — w, the above equations can be 
rewritten as 

D'a = {-1 + W + Xir)3nDd - 97^2(1 - w){l + h2f±^)Ei _ kUd + 9n{l - w)i\ir + A2 + 1 + w)$ 

I + w k 

+?,{Xir + \2)<^' - m7i{\ir + \2) - iUXirD^ , (14) 

U'd = 2{\ + -^(\ir + \2)\nUd + kDd-ik{\ir + X2 + l + iv)<^+{l + w)k^ . (15) 

By using the gauge-invariant quantity ( ~ (f) ~ Tidr and letting Q^n ~ Cd — C, ^c get the adiabatic initial condition, 

Dm ^ Dd , . 

1 - Ai - \2/r l + w + Air + A2 

The curvature perturbation relates to density contrast by [30] 

Al:Ga^Y.P^{Dl - P^U./P.J^ + m)k} 

k^-4TrGa^J2p'jH ' ^ ' 

With the help of these equations we can compute the curvature perturbation $ based on the CMBEASY code. 
We first consider the interaction between the dark sectors in proportional to the energy density of DM (A2 = 0) and 
in our calculation we keep the DE EoS w 7^ — 1. For constant w > — 1, we show the numerical results for the ratio 
r = Pm/pd in Fig.l. We observe that Air exhibits a scaling behavior, which keeps constant both at early and present 
times of the universe. This behavior is not changed when we turn on A2. Analytically, this can be understood by 
inserting the continuity equations 

p'„^ + 3Hpm = 3H{XlPm + X2Pd) 
p'd + 3Hpd{l + w) = -3H{Xipm + \2Pd) (18) 



^" r^ . (19) 
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Solving the corresponding quadratic equation, we get 



(rAi)i == -~-{w + Ai + A2) + -Jw"^ + 2u;A2 + 2wAi + A^ - 2A1A2 + A^ 
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(rAi)2 = --{w + Ai + A2) - -^w'^ + 2wA2 + 2w;Ai + A^ - 2A1A2 + A? . (20) 

This implies 

(rAi)i « — (W + A2) 

(rAi)2 « -^ - (21) 

w + A2 

for Ai ^ A2 < —w. These two roots of Air are constant in the very early time and current time of the universe, 
respectively. 

The scaling behavior of Air influences the curvature perturbation <I>. Numerically, we see from Fig.l that when 
w > —1 and Ai 7^ 0, $ blows up, which agrees with the result obtained in [27]. We find that this blow-up starts at 
earlier time when w approaches —1 from above and it happens regardless of the value of A2. 

The reason for the blow up is the fact that the expression of r is non perturbative in Ai, being proportional to Aj^ 
at very times, when we have to consider the bcgining of the CMB computation. 
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Figure 1: The upper two figures show the scahng behavior of Air. The lower two show the behavior of the perturbation. 



However, when we consider the dark sectors' interaction as being proportional to the energy density of DE and 
examine the case that the constant EoS is a httle bigger than —1: the blow-up disappears. StabiUty is also found 
when we extend our discussion to the constant EoS w < ^1. 

Numerically, we find that the first two terms on the RHS of eqs. (14) and (15) contribute more than other terms 
to the divergence. Using ^i and ^2 to represent the first two terms of eq. (14), we can approximately write 



where 



6 



Dd-^i+ 6 



-1 + w + \ir)3nDd , 

-9H^(l-.)(l + ^f±^)^ 
1 + w k 



(22) 



(23) 



When A2 = 0, Ai 7^ and — 1 < w < 0, Air sa — w, we have ^2 one order larger than ^1 and f 1 + ^2 > 0, which causes 
the vast increase in Dd as shown in Fig. 2a. However, when Ai ^ in the case w < —1 and Ai = 0, A2 ^ no matter 
whether w < —1 or w > —1, ^2 and ^1 are of the same order as shown in Fig. 2b, c and ^1+^2 < 0, which makes Dd 
to decrease with time. Therefore, the blow up is avoided. 

In order to further explain the reason for the blow-up wc provide an analytical analysis. Keeping the leading terms, 
we have the approximate equations 



D' 



U'd 



{-1 + W + Xir)3HDd - m^{l -w)il + bl±2^)l!± 

1 -f 1/7 k 



1 



1 



-(Air + A2) 



nUd + kDd. 



(24) 



Considering the case that the interaction between dark sectors is proportional to the energy density of DM (Ai 7^ 
0, A2 = 0) and noting that Air ~ — tu, we can simplify the above equations to 



Dl, 



U'd 



-3HDd - 9H- 

1 + W 



l-wUd 
1 + w k 

~kDd . 



A second order differential equation for Dd is 



D'^ « (2^ - \^w) D'd + m' n')Dd 



(25) 



(26) 
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Figure 2: Comparisons of ^i and ^2 



In the radiation dominated period, we have 7i ~ i,7i' ~ ^7^^ TT ^ ^t ^^"^ ^^' (^®) ''^'^ ^^ approximated as 



l + 3w;Z)^ 6 



(27) 



whose solution is 



Dd « Cir'-i + C2T 



(28) 
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where 



Figure 3: Behavior of the indices n and r2 in terms of w. 
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(29) 



The result, eq(29), has also been given in [27], see their equations (84) and (85) by setting a = 0, putting the obvious 
+/— in front of the square root and remembering D^, ~ ip. In [27], n corresponds to ip while r^ corresponds to 
D. In fig. 3, we see that when —1 < w < 0, both ri and r2 are positive, which means that the perturbation in D^ 
grows. However, when w < —1, both ri and r2 are negative; this results in the decay of the perturbation in Dd- No 
divergence occurs, regardless of the value of Ai. 

These results tell us that for a constant DE EoS w > —1, a, coupling between DE and DM in proportional to 
Pm(Ai 7^ 0) will lead to a violent divergence in the curvature perturbation. However, this divergence is absent for 
w < -1. 

Considering the case that the interaction between dark sectors is proportional to the energy density of DE, namely, 
Ai = 0, A2 ^ 0, eq(24) reduces to 



D'd w i-l + w)3nDd-9n^{l~w)[l + 
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(30) 



We can rewrite the second order differential equation for Dd in the form 
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(31) 



In the radiation dominated era, the above equation becomes 



D'^ = -3 + 3w 
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(32) 



Introducing the auxiliary quantities 
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Figure 4: Behavior of A and of F. 



we find that, when A > 0, 



where 



On the other hand, when A < 0, 
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In fig. 4 we see that A can be positive only in the vicinity of w = —1. When A2 is small, the range for positive 
A is small, w = — 1 is the central singularity, since it will lead to the divergence in ri and cause the blow-up in the 
density perturbation eq(34). When w > —1 and A > 0, the blow-up in the density perturbation can also happen 
since r/2(l + w) is positive as well. But when w grows further above —1, A will become negative and so does 
r/2(l -I- w), which will lead to the convergent result of eq(36). When w < —1, r/2(l -I- w) is always negative, the 
density perturbation will decay even when w is close to —1 from below and A is small and positive. 

The physical origin of such a behaviour can be traced to eq (8). When Ai = 0, the dark energy sound speed 
depends only on dark energy parameters, contrary to what happens when Ai ^ 0. In this latter case, the coupling 
introduces a dependence of the pressure perturbation on the dark matter energy density. In the latter case, at early 
times pm >> Pd and the non-adiabatic pressure perturbation diverges at super horizon scales, driving the instability. 
In our case, the effect is less acute and the system of coupled differential equations describing the evolution is better 
behaved. 

These results show that when the interaction between dark sectors is proportional to the energy density of DE(A2 7^ 
0), the blow-up in the perturbation will not happen for constant EoS w < — 1. For w > —1, when the coupling is 
small, the blow-up can also be avoided in the observational range of the EoS. However, there is a possibility for the 
divergence to happen when the interaction is large in the obscrvationally allowed w > — 1 range. 

In summary, we have reexamined the cosmological perturbations when DE and DM interact with each other. We 
have specialized the interaction to be a linear combination of DE and DM energy densities, namely XiPm + \2Pd- We 
found that for constant DE EoS w > —1 and nonzero Ai the instability occurs in agreement with the results of Ref. 
[27]. However when w > ~1 and the interaction is just proportional to the energy density of DE(Ai = 0, A2 7^ 0), 
the perturbation is stable for small A2 when w is within observational range. For phantom DE case with constant 
w < — 1, the perturbation is stable regardless of the value of the coupling. This result was also evidently shown in 



[27]. It would be interesting to extend this study to other interaction forms. Moreover, it would be of great interest 
to confront the stable DE and DM interaction model to observations, such as CMB angular power and large scale 
structure etc. Works in these directions are in progress. 
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